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Abstract. - The properties of geodesies flow are studied in a Friedmann- Robertson- Walker metric 
perturbed due to the inhomogeneities of matter. The basic, averaged Jacobi equation is derived, 
which reveals that the low density regions (voids) are able to induce hyperbolicity, even if the global 
curvature of the Universe is zero or slightly positive. It is shown that the energy independence 
is a characteristic property of these geometric effects. The importance of these conclusions is 
determined by the temperature independent ellipticity of excursion sets and regions of different 
randomness found in Kolmogorov CMB maps. 
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Introduction. — The increase in the accuracy of Cos- 
mic Microwave Background (CMB) radiation measure- 
ments has not only largely revealed the set of the basic 
cosmological parameters [1], but has also stipulated the 
importance of the search of non-gaussian effects. The list 
of already reported non-gaussianitics includes the low mul- 
tipolc alignment anomaly, the Cold Spot, etc (e.g. [2] 
and references therein). The CMB distortions due to the 
cumulative contribution of inhomogeneous matter distri- 
bution in the Universe are among the discussed effects, 
which on the one hand, can influence the primordial non- 
gaussianities, on the other hand, will enable the link with 
the observed filaments (e.g. [3]). 

Another empirical fact for potential interest here is the 
temperature independent ellipticity of anisotropics found 
in CMB sky maps [4]. For example, it was noticed [5] 
that, if the nature of the dark energy is due to the matter 
inhomogeneities, then it might be linked with those ellip- 
ticities in the CMB maps. For further consideration of the 
backreaction of the inhomogeneities see [6]. 

Here we consider within a geometric framework the 
CMB distortions both due to the matter distribution and 
the overall geometry. This approach is based on the con- 
cept of the instability of geodesic flows and involves the 
Hamiltonian formalism, known by its efficiency in the 
problems of instability of dynamical systems [7]. 

It is shown that the temperature independence of the 
distortions of CMB excursion sets is a generic feature of 



the geometric nature of the effect, as studied initially in [8]. 
Namely, whatever is the agent causing the hyperbolicity 
itself, one has an empirical indication of a geometric effect, 
thus supporting such a nature for temperature threshold 
independence of ellipticitics [4] . 

Then, we show that the low density regions, voids, can 
induce hyperbolicity even in conditions of globally flat or 
positively curved Universe. 

To have possibly self-consistent account, we start with 
a brief account of the reduction of the null geodesies from 
(3+1) to 3-spaces, then move to the averaged Jacobi equa- 
tion and the conditions for the instability of the flows. 

Anosov flows. — The basic steps can be separated 
as: (a) to find out the conditions when the geodesies in 
(3+l)D can be projected as geodesies in 3D Riemannian 
manifold, (b) the Hamiltonian representation of the flows, 
(c) the time correlations functions of Anosov flows. 

Let M be a smooth d dimensional manifold with a 
smooth Riemannian metric g, t : TM — > M be the tan- 
gent bundle over M, and K : TE — > E the linear connec- 
tion on a vector bundle 7r : E — > M over M, such that 
for each m € M and X € E m = 7r _1 (m) the restriction 
K : TxE — ► E m is linear and surjective. Then, the follow- 
ing diagram is commutative 



TE 

TE I 

E 



E 
M. 
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For X G E m we define 

Ver(X) = kcr(T7r : T X E T m M), 

and 

ffor(X) = kcr(K x : T X E -» £ m ), 

to be the vertical and horizontal subspace of TxE, respec- 
tively. Therefore, 

T X E = Hor(X) © Ver(X). 

The covariant derivative V of a differentiable section £ : 
M — > _E is defined by V£ = K oT^. For any smooth curve 
c : -R — ► M : t \— > c(i) with tangent vector c = Tc.^ = u 
and tensor field £ along the curve c we have 

The isomorphism between T X TM and T ToX M ®T ToX M 
can be stated as 

TtxK : T X TM T toX M®T toX M : Y ^ (Tr.Y, ifJ). 

It is well known that on a smooth Ricmannian mani- 
fold there exists exactly one torsion-free Riemannian con- 
nection, the Levi-Civita connection. Thus, for the Levi- 
Civita connection the Ricci identity holds \7[g(X, Y)] = 
g(VX,Y) + g(X,'VY), and it is torsion-free, i.e. for any 
vector fields X and Y we have V X Y - Vyl = [X,Y]. 
The Riemannian tensor is given in terms of a covariant 
derivative by the following formula 

Riem(X, Y)Z = ([Vx, Vy] - V [JC ,y])Z. 

The linear transformation ?R. X : Y — > Riem(Y, X)X is 
known as the curvature transformation. 

Let 5* be a smooth vector field on the tangent bundle 
TM of a smooth Riemannian manifold M 

S : TM -> TTM : u h-> («, F) G Hor(r o it) © Ver(r o u). 

We denote by u : R — > TM : i i-» = /*ito an integral 
curve of the vector field 5 passing through the initial point 
uq G TM i.e. u — S(u), and u(0) = uq. The group of 
diffeomorphisms {/'}, t G R 

f : TM -> TM : u ^ u(t), 

is a one-parametric dynamical system or flow of the vector 
field S on TM. The curve u(t) is the flow line or integral 
curve starting at u(0) = uq, c = r o u is a trajectory. 

If u(t) is an integral curve of the vector field S starting 
at uq, then c = r o u is a solution of the following equation 

V £ c = F, 

determined by c(0) — tuq and c(0) = uq. In addition, if 
V £ c = F, then c(i) = /*c(0). 



We will discuss vector fields S, such that i 7, = 0. In 
order to study chaotic properties of dynamical systems we 
need to investigate solutions of the equation for invariant 
vector fields Z along the curve c, i.e. Z(t) = Tf t Z(0). 
They are in 1 : 1 correspondence with the solutions of the 
Jacobi equation 

V 2 u z + $t u (z) = 0. (1) 

The correspondence is given by 

z(t) <-> Z(t) = (z(t),V u z(t)) G Hor(u(t)) ®Ver(u(t)). 

It can be proved that the dynamical system defined by 
a vector field S(u) = (u,0) is a geodesic flow with La- 
grangian L(x, u) = \g{u, u) 

c = u, 

< V u u = 0, (2) 

y 2 u z + $t u (z) = o, 

and g(u, u) = 1. The Jacobi field z (solution of the Jacobi 
equation) is said to be orthogonal to it, if g(u, z) = 0. Any 
Jacobi field has the following form z = n + const ■ u, where 
n is orthogonal to u, i.e. g(u, n) = = g(u,\7 u n). Thus, 
n is orthogonal to u for all t, if is orthogonal at t = 0. It 
is a very important property of the geodesic flow, for it 
means that equation for n is the same as for the z. 

One can derive an equation for the length of n from the 
Jacobi equation 

V 2 u n + n u (n) =0. 

Let n = in, where g(n,n) = £ 2 , g{h,h) = 1, g(h,u) = 0, 
then we get 

i + [K(u,h)-g{V u n,V u n)]i = 0, (3) 

where the sectional curvature has the following form 

K(u,h) = g(^R u (n),n) = g^R.(h,u)u,h). 

If M is compact and K(X,Y) < K < for all or- 
thonormal vectors X and Y at any point of M, then the 
geodesic flow is an Anosov system [9]. This condition is 
rarely fulfilled for physical systems. Nevertheless, consid- 
ering below the perturbed geodesic flow and to outline 
the link with basics of Anosov systems, we will refer as 
of Anosov instability, even though the inequalities are ful- 
filled for averaged quantities only. eq. (3) is usually very 
difficult to solve, so one is forced to find a simpler equa- 
tion, the solutions of which approximate the original one. 
Substituting 

1 / «®u\ 

d- 1 v g{u,u)) 

for n®n in K(u, h) we get the following averaged equation 
£+ Ric{u,u) £ = 0, (4) 



Flat Universe with Hyperbolic Voids 



where we have dropped the non-positive term 
— ff(V M n, V u n) < from eq. (3). One may average 
this again by replacing u <g> u with ^9 to get 



d(d-l) 



I? 







(5) 



(R is the scalar curvature). These operations are justi- 
fied by the search of averaged properties of geodesies in 
Robertson- Walker homogeneous and isotropic metric, and 
the fact that the solutions of the averaged equation are 
typically greater those of the original equation eq. (3) (cf. 
[7]). 

Reduction of (3+l)D Null Geodesic Flow to 3D 
Geodesic Flow. - Let V = M x R 1 be a (3+l)D 
smooth manifold with the following metric 



ds 1 



g flu dx fJ 'dx 1 ' 
a 2 ( V ) e 2 ^ 



-drf 



(6) 



where 



da 2 



j mn {x)dx m dx n 



is any smooth metric on M; \i, v = 0, 1, 2, 3; m, n = 1, 2, 3. 

As is well known the geodesic flow on V can be described 
by the following Hamiltonian 



H(p, x) 



\g^ u (x)p n p v 



\a- 2 



(v) e 



-2i/> 



-pl + e 2 f^j mn p m p„ 



and Hamiltonian equations 

dx>* _ dH 
ds dp^ ' 



i± 

ds 



dH 



Moreover, for the null geodesies we have H(p, x) = 0. 

We are going to prove that projections of the geodesies 
from V into M will be geodesies on M having the metric 
7 = e _2 ^7. It is easy to see that for the null geodesies 



dH = m = _ 2 & n 

dx n drj a 



0. 



therefore 



and 



ds 



= 0, 



dr\ 
ds 

thus, po is a constant (E 
nates we obtain 

dx m 



-a- 2 e^p , 



Po)- For the "space" coordi- 
dH 



where 



H 



ds 
ds 



\a- 2 e~ 2 ^ 



dp,, 



PrnPr 



-\a- 2 e~ 2 ^] 



1 „2f ran 



dH 

dx m 



If we define dr = a 2 e 2 ^ds, then we get 
dx m dH dp m dH 



dr 



dp-n 



dr 



dx m ' 



with H = \pq. Thus, the (3 + 1)D null geodesic flow 
(V, g, H, ds) is reduced to a 3D geodesic flow (M, 7, H, dr). 
This generalises results of [11]. 

Energy dependence of the mixing. It is obvious 
that 



a -2 e -W = <i l = dr ldT 



= [-a- 2 e- 2 ^p ] 



dr 
dr/ ' 



ds ds drj 

hence dr = —drj/po- For H = H/p 2 , and f = p^r we have 
dx m dH dp m dH 



df 



dp n 



dr 



dx r 



with 



where 



f-J — I n mn n n — I 

±± o.y PraFn o' 



9mn — Po 'Jran • 



From Ergodic theory (e.g. [10]) we know that the rate of 
the decay of correlations (also called the rate of mixing) 
is proportional to the e _XT , where Lyapunov exponent 



X~ VR(§) = \po\- 1 VW) and 



\po\f]- Hence, 



XT - 7] 



VW)- 



Therefore, the mixing rate does not depend on the en- 
ergy of photons (E = — po). Thus, the relevant observable 
effects, such as the ellipticity of the excursion sets [8], have 
to be independent on photon's energy, temperature. 

Instability in perturbed FRW Universe. Let 

V = M x R 1 be a (3+l)-D smooth manifold with a per- 
turbed Robertson- Walker metric of the form 

ds 2 = g^ydx 11 dx v 

= -{l + 2cj))dt 2 + (l~2(j))a 2 {t)do- 2 , (7) 



where 

Iran = (l + \ [^ ) 2 + ( X 2 ) 2 + (z 3 ) 2 ] ) ~* S mn 

is the metric of a 3-sphere (k = 1), a 3-hyperboloid (k = 
— 1), a flat 3-space (k = 0), respectively. As in [12] we 
assume that 

\<t>\ « 1, 

and time derivatives of </> are much smaller than spacial 
derivatives, i.e. 



and 



2<i 7 PraPa 



||V0|| 2 « |A0|, 
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where 



where 



dx m dx r - 

If we substitute ip = 4> and / = 20 and take into account 
that e 2 * ~ 1+20 and e~ 2 ^ ~ 1-20, then the metric cq. (6) 
will have the same form as eq. (7) for d = 3, therefore we 
can use the above results. 

To estimate the Lyapunov exponents one has to inves- 
tigate the Riemann scalar curvature for the reduced space 
metric 7 = e _4< ^7 (see eq. (5)). It can be easily shown 
that 



where [12] 
Therefore, 

Then, from 



R = 6e^ [fc-f(A0+||V0|| 2 )] 
-A0 = 47ra 2 <5p. 



-~fc--A0 = /c+ — a 2 *p. 



3 ' 



k_ A 87rp a| 
a 2 3 3 a 3 ' 



ao l + z a 



to 



dr 1 r 



t a(r) a J Q H(g) 



one has 



1 = n k + n A + Q m , 

H 2 = Hi [fi A + [1 - fi A + O m z] (1 + z) 2 ] 



where 



/C A 87T/9Q 



a 2 T 3 T 3 ' 



3-ffg 



a n Ii n orin on. 



Thus, 



f/(*o) - »?(*) = ( a o#o) 1 A(z, f2 A , fi m ), 



where 
X(z, 51a, n m ) 



y/n A + [1 - r» A + siu] (i + 



(8) 



(9) 



<5o 



A) ' 



Finally, the eq. (5) can be rewritten in the following form 



d\ 2 



+ r£ = 0, 



where 



r = + 25 Q r 



(10) 



(11) 



By definition Sq > — 1, therefore, r > — fi/. — 2fi m . 

These equations are the basic ones describing the mix- 
ing properties of null geodesies, i.e. of the propagation of 
CMB both due the global geometry and local perturba- 
tions of metric (lensing) . 

For example, when 5q = 0, then r = — Ofc, as investi- 
gated in [8]. While in the inhomogencous Universe, the 
resulting effect can be different. 

Indeed, the underdense regions with 5 < 0, the voids, 
contribute to the hyperbolicity. Importantly, this conclu- 
sion does not require compactness of M. 

Conclusions. — The geometric approach of theory of 
dynamical systems is used for the description of the insta- 
bility (mixing) of null geodesies peculiar to Anosov flows, 
both due to global curvature and locally perturbed metric 
in the context of the properties of CMB. The correspond- 
ing equation, of averaged Jacobi form, and the Lyapunov 
exponents are derived, from which it follows that the low 
density regions, i.e. the voids, would contribute to hy- 
perbolicity, even if the Universe is globally flat or slightly 
positively curved. 

The available observational parameters of the voids (i.e. 
of 30 Mpc size) obtained from the large scale galaxy red- 
shift surveys seem to support their role of hyperbolic lenses 
[13], namely, due to the cumulative effect of such voids 
(parametric resonance), and not due to a single void of the 
horizon scale. The hyperbolicity of the voids and hence 
the loss of correlations in the propagating photon beams 
can be measured by the assigned degree of randomness in 
the CMB maps [14]. 

It is also shown that the peculiar property of the geo- 
metric mixing is its independence on the energy of pho- 
tons. This fact is in agreement with the temperature inde- 
pendent ellipticity of anisotropies found in the CMB sky 
maps, first of COBE, then of Boomerang and WMAP [4], 
thus supporting the geometric nature of the effect. 

This approach can be efficient also at numerical simu- 
lations of CMB sky maps cross-correlated with the distri- 
butions of the filaments. 



And 



R 



{a Q H ) 2 {-Q k + 2J n m ) = (a iT ) 2 r, 
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TT\-A 



Flat Universe with Hyperbolic Voids 



REFERENCES 

[I] de Bernardis P. et al, Nature, 404 (2000) 955. Spergel 
D. et al, ApJS, 170 (2007) 377. Komatsu E. et al, 
arXiv:0803.0547. 

[2] Rakic A., Schwarz D.J., Phys.Rev., D75 (2007) 
103002. Copi C.J. et al, Phys.Rev., D75 (2007) 023507 
Gurzadyan V.G., Starobinsky A. A. et al, A & A, 490 
(2008) 929. 

[3] Lewis A., Challinor A., Phys.Rept, 429 (2006) 1. 
Slosar A., Hirata C. et al, arXiv:0805.3580. Granett 
B.R., Neyrinck M.C., Szapudi I., arXiv:0805.3695. 

[4] Gurzadyan V.G., Torres S., A & A, 321 (1997) 
19. Gurzadyan V.G., de Bernardis P. et al, 
Mod.Phys.Lett., A20 (2005) 893. Gurzadyan V.G., 
Bianco C.L. et al, Phys.Lett. A, 363 (2007) 121. 

[5] Wiltshire D.L., New J.Phys, 9 (2007) 377. Leith B.M., 
Cindy Ng S.C., Wiltshire D.L., ApJ, 672 (2008) L91. 

[6] Zalaletdinov R.M, Gen.Rev.Grav., 24 (1992) 1015. 
Buchert T., Gen.Rel.Grav., 40 (2008) 467. Buchert 
T., Carfora M., Class. Quant. Grav., 25 (2008) 195001. 
Larena J., Alimi J.-M. et al., arXiv:0808.1161. 

[7] Arnold V.I., Mathematical Methods of Classical Mechan- 
ics (Springer- Verlag) 1978. 

[8] Gurzadyan V.G., Kocharyan A. A., A & A, 260 (1992) 
14; Europhys. Lett., 22 (1993) 231. 

[9] Anosov D.V., Comm. Steklov Mathematical Inst., 90 
(1967) . 

[10] Handbook of Dynamical Systems, edited by Hasselblatt 
B., Katok A. (North Holland) 2002. 

[II] LOCKHART CM., MlSRA B. AND PRIGOGINE I., PhyS. 

Rev., D25 (1982) 921. 
[12] Holz D.E., Wald R.M., Phys. Rev., D58 (1998) 063501. 
[13] Gurzadyan V.G., Kocharyan A. A., A & A, 493 (2009) 

L61. 

[14] Gurzadyan V.G., Kocharyan A. A., A & A, 492 (2008) 
L33. Gurzadyan V.G., Allahverdyan A.E. et al, A & 
A, in press (2008) , arXiv:0811.2732. 



